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Abstract 

An experiment has been conducted to acquire data 
for the validation of computational fluid dynamics (CFD) 
codes used in the design of supersonic combustors. The 
primary measurement technique is coherent anti-Stokes 
Raman spectroscopy (CARS), although surface pressures 
and temperatures have also been acquired. Modern- 
design-of-experiment techniques have been used to 
maximize the quality of the data set (for the given level of 
effort) and minimize systematic errors. The combustor 
consists of a diverging duct with single downstream- 
angled wall injector. Nominal entrance Mach number is 2 
and enthalpy nominally corresponds to Mach 7 flight. 
Temperature maps are obtained at several planes in the 
flow for two cases: in one case the combustor is piloted by 
injecting fuel upstream o! the main injector, the second is 
not. Boundary conditions and uncertainties are adequately 
characterized. Accurate CFD calculation of the How will 
ultimately require accurate modeling of the chemical 
kinetics and turbulence-chemistry interactions as well as 
accurate modeling ol the turbulent mixing. 

Nomenclature 
C N2 concentration of nitrogen 
A\ thermal conductivity of wall (W/mK) 
ti number of samples 

p number of parameters 

q heal flux (W/m 2 ) 

t time (s) 

T temperature (K) 


v Vf j position coordinates in right-handed system, 
see Fig. 2 (m) ^ 

a x thermal dittusivity of wall (m /s) 

G standard deviation 

Introduction 

Computational fluid dynamics (CFD) codes arc 
extensively employed in the design of high-speed air 
breathing engines. CFD analysis based on the Reynolds 
averaged Navier-Stokes equations uses models lor the 
turbulent fluxes that employ many ad hoc assumptions and 
empirically determined coefficients. 1 ypically, these 
models cannot be applied with confidence to a class ol 
flow for which they have not been validated. This 
experiment is one of several adopted by a working group 
of the NATO Research and Technology Organization 
(RTO) as a test case for their CFD development and 
validation activity. Another is a study of supersonic 
coaxial jets without combustion'. Calculations of the 
coaxial jet and the present geometry using the SPARK 
code 2 , performed in conjunction with the RTO activity, 
have been presented’ . 

The experiment is designed to provide a relatively 
simple case for CFD codes involving supersonic in jection, 
mixing and combustion in a duct. The model geometry is 
simple and large regions of subsonic/recirculating flow 
are avoided. Care is taken to define the geometry and 
model inlet conditions. The enthalpy of the test gas (hot 
air •■simulant”) is nominally equivalent to Mach 7 flight. It 
was believed, on the basis of calculations performed (e.g.. 
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Rel. 2), that this would produce mixing-limited flow, that 
is to say, one for which chemical reaction to equilibrium 
proceeds at a much greater rate than mixing. It later 
proved that this was not the case. 

The primary experimental technique employed is 
coherent anti-Stokes Raman spectroscopy, known by 
its acronym CARS. An introduction to CARS is given 
by Eckbrelh 4 . and an application of CARS to 
supersonic combustors is given by Smith et al.\ The 
species probed is molecular nitrogen and the quantity 
measured is temperature, although. If assumptions are 
made as to the relative proportions of the other major 
chemical species, composition can also be estimated. 
Composition data are not reported herein. CARS, like 
other optical techniques, does not significantly alter the 
flow being studied. Intrusive probes, such as Pitot, total 
temperature, hot-wire, etc., are not used due to access 
difficulty and high heat flux in the combustor, and 
because they may alter the flow. Also, CARS has 
several advantages over other optical methods. It is 
relatively mature and well understood. Signal levels are 
relatively high. The signal is in the form of a coherent 
(laser) beam, and can be collected through relatively small 
windows. Consequently, incoherent (non-CARS) 
interference can be rejected by spatial filtering. 

Application of a complicated technique like CARS 
in high-speed engine environments is not routine. Since it 
is a pointwise (rather than planar) technique, building a 
“picture" of the internal temperature field of The 
combustor requires hundreds of facility runs, which is 
expensive. Thus, modern-design-of-experiments (MDOE) 
techniques are used to minimize the quantity of data 
required to meet the goals of this work. Due to the 
complexity of the experiment and the fact that, at the 
outset, the effects of important variables (including 
facility operation, model and instrumentation variables, 
etc.) are not fully understood, there are likely to be 
significant uncontrolled variables. MDOE techniques are 
used to minimize systematic errors associated with these. 

Flow Facility 

The experiment is conducted in NASA Langley's 
Direct-Connect Supersonic Combustion Test Facility 
(DCSCTF) . Hot air simulant, known as vitiated air, is 
produced in the “healer", shown in Fig. I . Oxygen and air 
are premixed and then hydrogen is burned in the oxygen- 
enriched air. Flow rates are selected so that the mass 
fraction of oxygen in the vitiated air is the same as that of 
standard air. The high pressure, vitiated air is accelerated 
through a water-cooled convergent-divergent nozzle 
before entering the test model. The facility is designed to 
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test the combustor of a supersonic combustion ramjet 
(SCRAMJET) engine by directly connecting the facility 
nozzle exit to the entrance of the combustor. This facility 
is typically used in fundamental combustor research. 



To model 


Heater 

liner Cooling 
water in 


Figure 1 DCSCTF heater and nozzle* 

The test conditions are nominally representative of 
Mach 7 flight. Gas flow rates to the heater are: 
0.9 1 510.008 kg/s air, 0.028410.0006 kg/s hydrogen, and 
0.30010.005 kg/s oxygen. The heater stagnation pressure 
is 0.76510.008 MPa. All uncertainties presented in this 
paper are based on the 95% probability limits ( 1 .98 times 
the standard deviation). The above uncertainties are due to 
the random run-to-run variations and do not include a 
13% uncertainty in the mass flow rate measurements. 

Heater and nozzle exit conditions are estimated from 
the flow rates, heater pressure, and nozzle minimum and 
exit areas using one-dimensional ( I D) analysis 7 . The flow 
exiting the heater into the nozzle is assumed to be in 
thermodynamic equilibrium, but has unknown enthalpy 
due to heat lost to the structure and cooling water. The 
enthalpy is guessed and the area of the sonic throat 
computed by 1 D analysis assuming isentropic flow in the 
nozzle. The enthalpy is then iterated until the computed 
area at the sonic throat equals the geometrical minimum 
area of the nozzle. Nozzle exit conditions are computed 
from the geometrical exit area. The composition at the 
nozzle throat and exit could be evaluated at frozen (at 
heater values) or equilibrium conditions. All significant 
minor species are included. Calculations assuming 
equilibrium and frozen composition differ in minor 
speaes concentration, but not significantly in major 
species, temperature or pressure. The nominal calculated 
conditions, and uncertainties due to mass flow rate 
measurement error and run-to-run variations in heater 
conditions are: heater stagnation temperature 1 827175 K. 
exit temperature 1 187160 K, exit pressure 10011.5 kPa 

A . , . exit Mach number 1.98910.005. Errors arising in the 
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calculation due to the assumption of 1 D flow (the effects 
of non-uniform composition, boundary layers, etc.) ate not 

considered. . 

A study of the flow quality at the exit of the facility 
nozzle was previously conducted 8 . A Pitot probe rake was 
employed to map the exit Pitot pressure and additionally 
the flowfield at the exit of the nozzle was visualized. 
Silane (SiH 4 ) was added to the heater hydrogen and 
burned to form silica particles in the heater. The particles 
were illuminated by a pulsed laser-sheet and imaged with 
a CCD camera. Results were compared to CFD 
calculations of the nozzle flow. The flow at the nozzle exit 
was not completely ID. but the computed Pitot pressure 
distribution agreed well with measurement. The flow 
appeared well mixed. 

The test model is shown in Fig. 2: flow direction is 
from left to right. The model consists ot two main sections 
of duct: the copper upstream section and the carbon steel 
downstream section. Stainless steel flanges and carbon 
caskets separate the sections from each other and the 
nozzle. Proceeding from left to right, there is a constant 
area segment, a small outward step at the top wall, a 
second short constant area segment followed by a constant 
3° divergence of the top wall. The span is constant at 
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Figure 2 Test model: (a) nozzle, copper and steel duct sections, (b) detail 
in vicinity of fuel injector and pilots. 
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87.88 mm. Five small pilot fuel injector holes are located 
ahead of the step, and the main fuel injector is located just 
downstream of the start of the 3° divergence. The 
injection angle is 30° to the opposite wall. The injector 
nozzle is designed by the method of characteristics to 
produce Mach~2.5, ID flow at the injector exit. Hydrogen 
injection is provided at a pressure of 2.1210.07 MPa . 
temperature of 302±4 K. and equivalence ratio of 
0 99±0.04. On some inns, additional hydrogen injection is 
provided by the 5 pilot injectors at the same nominal 
temperature and a total equivalence ratio ot 0. 1 48 ±0.008. 
The pilots are turned on and off at the same time as the 
main fuel injector. 

The duct is uncooled; however, the wall thickness ot 
the copper duct is greater than 32 mm and the carbon steel 
duct is 19 mm. Thus, given the good thermal conductivity 
of these materials, it is possible to operate the facility with 
the model fueled for run times in excess of 20 s (and 
unfueled for much greater times). With atmospheric 
temperature air flowing in the model between runs, runs 
could be repeated every 10-15 minutes. 

The model is equipped with 7 slots to allow the 
CARS beams to penetrate the duct, of which slots 1 . 3. 5. 

6, 7, depicted in Fig. 2(a). are used in 
this study. The slots axe in paiis, one 
on each side of the duct. 4.8 mm wide, 
extending the full height ot the duct. 
When not in use the slots are plugged 
flush to the wall. Windows covering 
the slots are mounted at the end ot 
short rectangular tubes, at the Brewster 
angle to minimize reflections (Fig. 3). 
The window tubes are ventilated with a 
constant flow of electrically heated (~ 
400 K) dry air to prevent condensation 
of water on the windows. The CARS 
beams focus can thus be translated the 
full span and height of the duct without 
damaging the windows. 

The model is instrumented with 
both pressure taps and wall 
temperature probes. The details of this 
instrumentation are different in the 
upstream and downstream duct 
sections. Thirty -five static pressure 
taps are located in the copper duct, 
consisting ot 0.80 mm diametet 
square-edged holes. Taps are located 
on the bottom wall, at the centerline, 
and on the top wall at ; = -36.3 mm (: 
is measured from the horizontal center 
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line). Forty-nine static pressure taps are located in the 
steel duct, consisting of 1.6 mm diameter square-edged 
holes. Taps are located at top and bottom wall centerlines, 
and sidewall midpoints. Pressures are measured with an 
accuracy of ±0.6 kPa by a Pressure Systems Inc. 
electronically scanned pressure-measuring system. 

Six spring-loaded, bayonet, ribbon, K-type 
thermocouples from Nanmac Corp. are located in 6.35 
mm diameter blind holes in the copper block. The bottom 
of these holes is square and the wall thickness (to the duct 
interior surface) is 2.8 mm. The thermocouple junction is 
located at the bottom of the hole at the center. This 
diameter and wall thickness are chosen, based on a 
multidimensional heat transfer analysis, so that the 
thermocouple would measure the unperturbed duct surface 
temperature. While the temperature at the bottom of the 
hole is less than local duct surface temperature, the effect 
of the hole, which alters the conduction path, is to raise 
the local duct surface temperature, and these effects 
offset. Six "eroding" K-type thermocouples from Nanmac 
Corp., 7.9 mm diameter and incorporating a carbon-steel 
sheath, are located in the carbon-steel duct. The 
thermocouple junction for these probes is located flush to 
the duct flow surface. Over time, the junction is removed 
by the hot flow’, but can be regenerated using sandpaper. 
Since the probe material is predominantly carbon-steel, 
the same as the duct, this probe also measures the 
unperturbed duct surface temperature. 

CARS Technique 

Optical System 

The CARS system uses an unseeded Spectra-Physics 
DCR-4 pulsed Nd:YAG laser frequency doubled to 532 
nm. The nominal power is 550 mJ per pulse in the green, 
repetition rate is 10 Hz, and nominal line width is less 
than I cm . Beams are horizontally or "p" polarized with 
respect to the plane of the optical system. A broadband 
dye laser is employed in the system. This laser consists of 
an oscillator cavity, formed between a total reflecting 
mirror and a 30% reflector, and a single amplifier stage. 
Two identical dye cells are used, one in the oscillator and 
one the amplifier, through which flows rhodamine 640 in 
methanol. After separating the 532 nm from residual 1064 
nm with dichroie mirrors, the 532 nm beam is split. One 
beam with 45% of the total power is used to pump the dye 
laser, split 22% to the oscillator and 78% to an amplifier. 
Pump beams are loosely focused on the dye cell, nearly 
parallel to the dye laser beam. Cell windows are at the 
Brewster angle to minimize reflections and ensure the dye 
laser is "p" polarized. The wavelength is centered between 
605 nm and 606 nm by adjusting the dye concentration. 
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matching the Raman shift of nitrogen. 

Figure 3 shows a schematic of the remaining 
components of the CARS system. The second 532 nm 
beam is split into two parallel, roughly equal intensity 
beams. The dye beam is expanded in a beam expander 
(B). Dye and 532 nm beams are combined at a dichroie 
mirror (C) and are relayed via a periscope to a spherical 
(focusing) lens (SL) of focal length 0.41 m. All three 
beams cross at their focal points: the focal point of the dye 
beam is made coincident with the focal point of the 532 
nm beams by adjusting the beam expander. Beams are 
phase matched in a vertical planar BOXCARS 
configuration 3 . 

At the lens, the full (to - 1 0% of peak) diameter of 
the 532 nm beams is - 8.5 mm, and the beams are 1 8 mm 
apart. The full diameter of the dye is - II mm. At the 
locus, the diameters are respectively - 0. 12 mm and ~ 
0.1 5 mm, yielding a geometrical beam intersection length 
of 5.4 mm. The length of the measurement volume is 
measured by translating the CARS measurement volume 
through a thin planar jet of nitrogen, surrounded by a 
coflowing jet of helium. The length over which CARS 
signal is recorded is — 4.5 mm and the lull width half 
maximum (FWHM) of the signal distribution is - 2.25 
mm. The measurement volume length is selected (by 
selecting the spacing of the beams on the lens) as a 
compromise between the desire to have high spatial 
resolution and high signal power. The following beam 
energy levels per pulse are obtained at the focusing lens; 
~85 mJ for each of the green, from 12 mJ to 24 mJ for the 
dye. 

The beams (including the CARS signal beam) are 
relayed via a second spherical (collimating) lens and a 
second periscope back to the optical bench. The 
overlapping CARS and 532 nm beams are separated in a 
splitter (S). The splitter consists of two 100 mm long 
dichroie mirrors that reflect -99.5% of the incident CARS 
signal while reflecting only -20% of the 532 nm, 
transmitting the rest. The CARS and 532 nm beams enter 
the splitter at 45° to the mirrors and undergo six back-and- 
forth reflections. The beam is then directed through 
additional filters as needed to reject residual 532 nm or 
for CARS signal attenuation. It enters a polarizer that 
allows only p-polarized light to pass (not a critical 
component since it is found that the signal is already well 
polarized). It is then focused by a pair of cylindrical 
lenses (CL) and enters a 1 m monochrometer with 1200 
groove/mm grating via an aperture (E). An EG&G PAR 
model 1420 intensified, linear, self-scanned silicon 
photodiode array detector (IPDA) is mounted at the exit 
plane of the detector. The detector consists of 1024 
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elements. 2.5 mm high by 25um wide, ot which the 
central 598 elements are used. The cylindrical lenses are 
adjusted to produce a horizontal focus at the detector 
(maximizing image sharpness) and a loose vertical focus 
(to ensure the detector provides a linear response to 
signal 1 *). An optical splitter (X) is located in front of the 
detector 10 . The splitter creates a secondary signal on the 
detector, identical to the primary but offset by 290 pixels 
and of 6.1 % the intensity. When the intensity ot the 
primary signal exceeds the dynamic range ot the detectoi. 
the secondary signal is used for analysis. 

The two top prisms of the periscope are mounted on 
stepping motor driven vertical translation stages. The two 
bottom prisms and the vertical translation stages are 
mounted on similar horizontal stages. By translating the 
vertical and/or horizontal stages in tandem (maintaining 
alignment) the measurement volume could be moved in 
the v and/or : direction. It is important that the beams are 
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Figure 3CARS system: beam combining optics, intersection and signal 
detection optics. 

parallel to the direction of motion of the translation stage. 

A reference, low-power helium-neon laser (He-Ne ) beam 
is directed through the periscopes, parallel to the direction 
of motion of both sets of stages, and through the center of 
the spherical lenses. The Nd:YAG and dye laser beams 
are periodically realigned by ensuring that they intetsect 
the He-Ne beam at their focus. 

CARS data acquisition is under the control of a 
personal computer (PC). The PC provides a continuous 
sequence of pulses to trigger the laser and synchronize the 
The I PDA and translation stage controllers are 
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initialized prior to data acquisition. When data is to be 
acquired, a control pulse is sent to the I PDA and 
translation stage controller, which then execute 
preprogrammed sequences. The IPDA scans at 30 Hz: 
scans 1. 4, 7 etc. are clearing scans and are not saved, 
scans 2. 5. 8, etc. are prescans, and scans 3, 6. 9, etc, ate 
data scans, coincident with the laser pulse. After 
completion of acquisition, resultant IPDA data is 
transferred to the PC. Two types of acquisition are 
employed. In the first, data is acquired at a single point in 
space. In the second, data is acquired while either the 
vertical or the horizontal stages are in constant velocity 
motion. 

CARS data are acquired in the supersonic combustor 
during multiple sets of test runs. During a set of runs 
(which might last as long as 5 hours), access to the model 
and optical system is prohibited for safety reasons. Test 
runs consisted of approximately 5 s during which the 
heater is operating but no fuel is 
injected in the model, followed by 
from II s to 20 s during which fuel is 
injected. CARS data is acquired over a 
period 2 s shorter than the period of 
fuel injection. Immediately after a run. 
10 s of data is acquired with the system 
operating as before, the dye laser beam 
blocked by a remotely operated flag. 
These ‘"background" scans measure 
non-CARS interferences such as 
scattered laser light. Dye laser spectra 
are acquired simultaneously with each 
data run by sampling with a fiber optic 
probe coupled to an Instaspec l A m 
spectrometer with a CCD detector, 
controlled by a separate PC. 

Just before and just after a set of 
test runs, “reference" CARS spectra 
are acquired in a low speed jet ot “PC- 
Duster" refrigerant gas. Since this jet 
contained no nitrogen, and is ot 
constant CARS susceptibility, the 
spectra reflect the spectral variation in dye laser power. 
Additionally, spectra are acquired with room temperature 
air flowing in the duct. 

The test facility poses problems in relation to the 
operation ot the CARS system. The effect ot acoustic and 
structural vibration is not significant. Howevei. laige 
temperature swings are common since the room in which 
the apparatus operated is continuously ventilated with 
external atmospheric air: additionally significant heat is 
radiated from the model. Thermal expansion ettects on the 
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mechanical components of the optical system are not 
sutficient to require remote adjustments (adequate signal 
levels are typically maintained for several hours). 
However, thermal expansions significantly change dye 
laser center wavelength and power. 

Data Reduction 

CARS data are analyzed on a separate 
workstation. Prescans are subtracted from data scans. 
Background scans (after subtraction of prescans) are 
averaged and subtracted from data scans. Both primary 
and secondary (produced by the splitter) CARS signal are 
contained within the data scan. If the primary is saturated, 
the secondary is selected for analysis. Data scans are 
divided by the reference spectrum to remove the effect of 
the dye laser spectral power distribution, and normalized 
to unit area (primary or secondary). Data are compared to 
a library of similarly normalized theoretical spectra to 
determine the temperature and nitrogen concentration. 
The pixel location of the stall of the theoretical spectra is 
allowed to float ±10 pixels from a nominal position in 0.3 
pixel increments. The combination of temperature, 
concentration, and pixel location that produces the least 
mean square deviation between theory and data is 
selected. 

Theoretical CARS spectra are generated using the 
program CARSFT 1 '. The combustion gases are assumed 
to be a mixture of nitrogen and non-resonant buffer gas. 
both having non-resonant susceptibility of 8.5x10 lfi 
cm' /erg. The nominal static pressure is assumed to be 1 
atmosphere. In fact, as may be seen in Fig. 7, the pressure 
at the CARS planes is as much as 20 kPa below and 35 
RPa above. The error introduced by this assumption is 
estimated by calculating theoretical spectra over a range 
ot pressures and temperatures and fitting to them 
assuming 1 atmosphere pressure. The difference between 
fit and true temperatures (fit minus true) is about 25 K per 
10 kPa difference from 1 atmosphere for Ts near 2000 K. 

I he Exponential Gap Model for collisional narrowing of 
the Raman line shape is used: a Voigt model was also 
tried and found to make little difference for the test 
conditions. A 532 nm laser line width of 1 cm ! is 
assumed: spectra are not sensitive to 532 nm laser line 
width tor the test conditions. An experimentally 
determined instrument probe function is used. 
Temperatures from 150 K to 3000 K in increments of 25 
K, and 29 non-uniformly spaced concentrations of 
nitrogen are used in the library. Typical spectra at several 
temperatures and concentrations are shown in Fig. 4 
giving an indication of the sensitivity of the spectra to 
these variables 
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The peak of the reference CARS spectrum shifts 
significantly during a set of runs. This shift can be 
attributed largely to a shift in the dye laser wavelength. 
However, there are other factors. Thermal effects on 
mechanical components of the optical system, beam 
steering by the hot gases in the duct, and even translation 
of the periscope stages all can change the position of the 
CARS beam focus at the entrance slit to the 
monochrometer, and hence the position of the spectra on 
the detector. Two techniques are employed to correct for 
this problem. The first is to allow the pixel location of the 
start of the theoretical CARS spectra to float in the fitting 
process, as described. The second is to estimate the 
reference spectrum from the run data itself. The technique 
is inspired by the observation that the directly measured 
dye-laser spectra could be accurately fit to a Gaussian 
function. 

The procedure is illustrated in Fig. 5. The first step 
is to analyze the data from a given run using the reference 
spectrum obtained with the refrigerant gas at the end of a 
set of runs. For a particular laser pulse, spectrum A is the 
best-fit theoretical CARS spectrum, B is the measured 
CARS spectrum (prior to normalization by the reference 
spectrum), and C, the ratio of B to A, is an estimate of the 
reference spectrum for this particular shot. All such 
estimates for a given run are averaged, and then fit to a 
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and the measured temperature compared to a 
calculation based on measured flow rates and 
equilibrium chemistry (including minor species). 
The mean temperatures for several 10 s (100 
laser pulse) runs are shown in Fig. 6. The 





spectrum. This procedure is typically iterated twice, with 
the latest best-fit theoretical spectrum used to determine 
an improved reference spectrum. 

The techniques used for acquisition and analysis of 
CARS data in the supersonic combustor are tested in a 
“Hencken” adiabatic, flat- flame burner burning hydrogen 
in air. Equivalence ratio (ratio of hydrogen rate to 
stoichiometric hydrogen rate for given air flow) is vaiied 


adiabatic flame 

The measurements agreed satisfactorily with 
calculation at an equivalence ratio less than or equal to 
one. However, for hydrogen rich flames, the measured 
temperature is. on average, nearly 150 K high. It is 
suspected that this is a problem with the way the CARS 
spectra ate modeled in the presence of excess hydrogen, 
rather than with the experiment. For example, the 
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hydrogen S( 1 0) rotational line is located at a Raman shift 
of 2267.5 cm 1 (compare Fig. 4) 12 . Hydrogen is not 
specifically included in calculating the spectral libraries. 
However, comparison of measured spectra with calculated 
spectra that include the line suggest that its presence is not 
alone sufficient to explain the discrepancy observed. This 
problem is being investigated, but for now it is accepted 
that measured temperature may be -150 K high in 
hydrogen rich regions. 

Modem Design of Experiments 

Randomization 

When repeated measurements or “replicates” made 
within a relatively short interval of time are more alike 
than otherwise identical replicates acquired over a longer 
time period, the observations are not independent. This 
condition is common in real experiments, and is caused by 
often unrecognized sources of systematic variation that 
persist over time, such as instrumentation drift, 
temperature effects, operator fatigue and learning effects. 
Sample means and variances estimated under such 
conditions are biased by generally non-reproducible 
functions of the systematic variation, which render the 
experimental results difficult to duplicate and of dubious 
value for confirming the predictions of computational 
models. Standard references on experiment design treat 
this phenomenon in detail. 13 ,4 ,5 

Fisher 16 recognized that independent observations 
occur less frequently in nature than generally assumed, 
and proposed the simple act of randomizing the set-point 
levels of independent variables as an effective way to 
ensure independence in case nature neglected to do so in 
some inopportune phase of an experiment. For example, if 
a horizontal scan of temperature is planned for several 
unilormly spaced locations across an operational 
supersonic combustor model, Fisher would recommend 
that the spatial coordinates be selected in random order 
rather than systematically. If the spatial coordinates are 
selected in a monotonically increasing sequence, starting 
on the left side of the supersonic combustor duct and 
proceeding systematically toward the right side, for 
example, wall-heating effects could cause measurements 
made on the right side of the duct to be biased high 
relative to measurements made on the left. This would 
impose an artificial asymmetry on the spatial dependence 
of temperature across the duct. However, if the spatial 
coordinates are selected in random order, some of the 
measurements made on the left would occur earlier and 
some would occur later, and similarly for the right side, 
ensuring that the effect of all time-varying systematic 
errors (not just wall effects) would then be independent of 
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spatial coordinates. This would defend against artificial 
asymmetries that might otherwise occur. For these 
reasons, spatial set-points are set in random order in this 
experiment. 

Response Surface Methods 

Response surface methods are applied in this study, 
in which the spatial dependence of temperature is 
represented in each of a number of planes in a supersonic 
combustor duct by a mathematical function of the spatial 
coordinates. This method has the potential for providing a 
very compact representation of the temperature 
distribution, by replacing a potentially large volume of 
individual measurements with a relatively small set of 
numerical coefficients that can be used to adequately 
predict the temperature for any combination of spatial 
coordinates within the range examined in the experiment. 
The method has the additional virtue that it provides 
estimates of the response variable (temperature) for 
combinations of the independent variables (coordinates) 
that are not physically set in the experiment. 

Because the response models are based on 
experimental data afflicted with uncertainty, predictions 
made with these models will likewise feature some degree 
of uncertainty. The specific uncertainty associated with 
any one model prediction depends on the variance in the 
response data used to develop the model, the number and 
location of points from which the model is fitted, the 
nature ol the model itself (number and form of individual 
model terms), and the specific coordinate location for 
which the prediction is to be made. However, Box and 
Draper show that the average variance across points used 
to fit the model is simply p<rhu where p is the number of 
parameters in the model, n is the number of points used to 
fit the model, and a is the standard deviation of those 
data. 

An inverse relationship between data volume and 
inference error risk characterizes all experimental 
investigations. In this case, for a given functional form of 
the model (which determines /?), and a given measurement 
environment (which determines cr), the average model 
prediction uncertainty can be driven to arbitrarily low 
levels by selecting a large enough volume of data, n. 
Conversely, by establishing in advance the levels of 
uncertainty that are acceptable, the volume of data 
necessary to achieve stated precision goals can be 
computed. This is useful in the design of an experiment, 
because it can reveal in advance when insufficient 
resources are available to acquire data in a volume 
necessary to achieve minimum precision requirements, 
and it can also prevent the acquisition of substantially 
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more data than necessary to ensure acceptable levels ot 
inference error risk. In the present study, for example, 
standard deviations in supersonic combustor temperature. 
a, as high as 300 K. are typical. (It should be emphasized 
that this high cr is due primarily to real turbulent 
fluctuations in the flow rather than instrument precision.) 
A typical response model might have 20 terms, p. To 
ensure an average standard error for the response model 
that is no greater than 25 K. say, a 2880-point sample of 
data would have to be acquired in each plane. This may be 
as much as an order of magnitude fewer points than would 
be acquired if response surface methods are not 
employed. 

Data volume estimates are increased when resolution 
requirements are added to the precision lequiiements (e.g. 
minimum temperature differences that must be resolved), 
but the same inverse relationship exists between data 
volume and inference error risk. Several excellent 
references provide additional details on response suiface 


Results 

Surface Pressure and Temperature 

Surface pressure and temperature data are piesented 
for two typical runs, one in which the pilot injectors aie 
operating and one in which they are not. 1 hese tuns aie 
chosen because the gas flow rates to the heatei, injectot 
and pilot and the heater pressure are all very close to the 
respective averages over the total set of tuns. 
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injection commences at time t| = 6.0 s (piloted) ot 6.4 s 
(unpiloted). Pressures vary widely in the upstream region 
due to the complex nature of the shock wave system 
created by the injectors and step. The pressure tor the 
piloted case is higher than the unpiloted case between the 
pilot (x = 0.074 m) and about x = 0.7 m. due to 
combustion of gas from the pilot and main fuel injector. In 
the unpiloted case, pressure generally falls moving 
downstream due to divergence of the duct, until 0.5 m 
where it rises rapidly, peaking at about 0.75 m. 
Presumably, there is minimal combustion upstream ot this 
region. Downstream of 0.75 m the pressure drops 
smoothly in both cases but is higher in the unpiloted case, 
despite the greater total injected fuel rate in the piloted 
case. Differences between the two cases suggest that 
significant combustion of the fuel does not take place 
upstream of 0.5 m in the unpiloted case, but that the 
combustion then proceeds to completion by 0.75 m. It is 
not expected that fuel and air are fully mixed at 0.75 m, so 
that further mixing and combustion is expected to occur 
downstream. Comparison between measurements at 10 s 
and at 22 s reveals only small differences. There is no 
suggestion that the combustion delay experienced by the 
unpTloted case is affected by the increase in surface 
temperature that occurs during the course of a tun. 
Surface pressure varies as much as 20 kPa oi moie 
between top. bottom, and sidewall taps (not shown) 
upstream of about 0.65 m. consistent with the effects of 
shock waves. Downstream of this point, there are no 
differences between the walls. 



Figure 7 Surface pressure distributions along 
centerline of bottom wall. 


Surface pressure distributions for the pressure taps at the 
bottom wall centerline, averaged over 1 s intervals, are 
shown in Fig. 7. Data are shown at 10 s and 22 s into the 
run: the heater is initiated at time / (1 = 1 s and fuel 



Figure 8 Wall temperature history, unpiloted, top 
wall: data and fits. 


Since the duct is uncooled, surface temperatures 
vary greatly during the course of the run. In the copper 
section, temperature is typically about ~ 360 K at the start 
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of fuel injection but rises to as high as ~ 610 K. In the 
carbon steel section, it typically is ~ 440 K at the start and 
as high as ~ 950 K at the end. These var iations in surface 


temperature are least-squares fit to the solution for wall 
temperature of a semi-infinite body at initially uniform 
temperature T.. subject to steps in surface heat flux from 
zero to </„ at time r„ and from q„ to q» + </, at t,: 
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The fit parameters are T„ </„, and q, and the tit is 
conducted out to t = 1 1 s. 
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Figure 9 Surface heat flux distributions: top and 
bottom walls, centerline 


Representative temperature histories and fits (top 
wall, unpiloted) are shown in Fig. 8. Fits often diverge 
from the data beyond 1 1 s. indicating that heat flux 
continues to change slowly during the run. The material 
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Large increases occur in the piloted case on the top wall, 
downstream of the main injector (.v = 0. 1 66 m). indicating 
the presence of combustion products near the wall. In the 
unpiloted case, a smaller heat flux rise occurs ahead of 0.5 
m (where the pressure rise starts), indicating either the 
start of combustion near the wall, or an increase in heat 
transfer coefficient. An increase in heat transfer 
coefficient would be expected to occur under the 
injection-induced streamwise vortex pair. 

CARS Temperatures 

Data is acquired over 201 facility runs over 10 test 
days. Except for one day. when laser beams clipped the 
edge of the duct window slots due to thermal expansion 
and movements of the duct, the vast majority of the data 
are found acceptable and analyzed. Data are acquired at 6 
oi 7 fixed points near the horizontal centerline. 6 s with 
only the heater operating followed by 18 s with both 
heater and fuel injected. Data are also acquired during 1 6 
s of horizontal motion of the translation stages at 5 mm/s. 
oi from 9 s to 1 8 s of 5 or 6 mm/s vertical motion, during 
which time fuel is continuously injected. Figure 10 is a 3- 
dimensional cutaway view of the duct with the locations 
of the data points overlaid. 

Ideally, each consecutive data point would be 
acquired at a randomly selected location and plane. Real- 
life constraints prevent this. Motion from one point to the 
next takes significant time: if it is required to move 
between randomly selected points during the course of a 
limited-time facility run, the number of data points 
acquired would be considerably reduced. Changing from 
one plane to the next required shutdown of the facility and 
hardware changes. Thus, all data at a given plane are 
acquired consecutively. Data are acquired either at the 
fixed points or during continuous constant velocity 
motion, as described. However, the order of this data 


property JaJ/k s is taken to be 36.7 kW K/mV' : for the 

copper duct and 1 2 kW K/mV ,: for the steel duct. These 
values have not been verified experimentally, so this 
analysis should not be relied upon except in a relative 
sense (i.e.. case to case, location to location within the 
copper duct, injection to no injection). Heat fluxes are 
presented in Fig. 9. both top and bottom walls, on the 
centerline. Heat flux varies in the range 1 .0 MW/nr to 0.3 
MW/m : before injection (</„). and in The range 0.7 MW/m’ 
to 1.8 MW/m- with injection (</,, + </,). Note that, before 
injection, heat flux at a given location varies 10% to 30% 
horn case to case, probably reflecting variation in the 
initial temperature of the wall between runs. With 
injection, heat flux at downstream locations shows a large 
increase relative to before 


acquisition and directions of the individual motions is 
randomized. 

The CARS temperature data are all fit to a cosine 
series bivariate function of order 5 at plane 1 and order 6 
at the other planes, with the number of fit parameters 
respectively 21 and 28. Fits are performed using 
commercial software" 1 . The number of data points 
acquired per plane is from 2000 to 4000 and the standard 
deviation of the data from the fit at the various planes 
ranges from 196 K to 304 K. Thus, the fitted functions 
tepresent an estimate of the mean temperature distribution 

with mean uncertainty, given by 1 .98 Oyfp/n , from 36 K 

to 59 K depending on the plane. It is important to point 
out that the uncertainty in the surface fits to the data is 
lower near the center of the measurement plane and higher 


injection due to combustion. 
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near the edge. This uncertainty does not include the effect 
of the previously noted error in fuel rich regions of the 
flow, or any other non-randomized error. 

Piloted k 



Figure 10 Cutaway views of data acquisition locations. 

After obtaining the estimated mean temperature 
surfaces, the square of the deviation between 
measurement and mean at each of the measurement 
locations is evaluated. These squared deviation points are 
similarly fitted to a cosine series bivariate function of 
order 5. Standard deviation of these fits ranges from 
78000 K 2 to 190000 K 2 . Thus, the fitted function 
represents an estimate of the mean squared deviation 
distribution with uncertainty from 15000 K to 32000 K . 
Or. taking typical values of uncertainty and mean square 
deviation then forming a percentage, uncertainty of order 
±25000/300 2 xl00 = ±28%. It should be emphasized that 
this large uncertainty is purely as a result of insufficient 
data points, i.e., as n gets large it will become small. An 
additional error in mean squared deviation arises from 
CARS instrument precision. Expressed as a standard 
deviation, it was given as - 100 K. Thus, the experimental 
fitted mean squared deviation surface is everywhere of 
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order 10000 K 2 high due to instrument precision. Given 
these high levels of uncertainty and error, the mean 
squared deviation surface fits may be only qualitatively 
useful. 

Figure 1 1 contains 3-dimensional cutaway views of 
the duct showing contour plots of the fitted temperature 
functions (left) and the square root of the fitted mean- 
square deviation functions, or rms (right). Flow direction 
is from top left to bottom right. Recall that the main fuel 
injector is on the top wall between planes 1 and 3, and 
that the pilot injectors are at the top wall upstream of 
plane 1. 

Consider first the fitted mean temperature surfaces. 

At plane 1 in the unpiloted case the temperature is fairly 
uniform, between 1030 K and 1250 K. The mean 
temperature for all the data points of this plane and case is 
1162 K. This mean compares favorably with the value 
computed assuming 1-D flow from the heater, which is 
1 1 87±60 K. The fitted temperature is also reasonably 
uniform in the piloted case, despite the injection of pilot 
fuel upstream of this plane, except that the temperature 
drops close to the top wall, where the (cold) pilot luel is 
injected. There is no indication of pilot fuel combustion in 
this plane. 

At plane 3, in the unpiloted case, there is a region of 
low temperature which is the injected fuel plume, with 
temperatures as low as ~ 250 K at the center. There is no 
evidence of combustion of the injected fuel at this plane. 
In the piloted case, there is a band of hot pilot fuel 
combustion products close to the top wall. The center of 
main fuel plume may be seen as a cool region with 
temperature as low as - 650 K, greater than in the piloted 
case, suggesting some combustion. Plane 5 is similar to 
plane 3. Temperatures have risen near the center of the 
fuel plume, to as low as ~ 550 K in the unpiloted case and 
~ 1250 K in the piloted. The height of the region of hot 
combustion products near the top wall in the piloted case 
has increased. At planes 6 and 7 in the unpiloted case, 
temperatures have risen abruptly as compared to plane 5, 
suggesting nearly complete combustion (i.e., combustion 
which, if locally fuel rich, has consumed all available 
oxygen, or if fuel lean, has consumed all available fuel). 
The hot region close to the top wall (temperatures as high 
as ~ 2300 K), and to a lesser extent that near the bottom, 
are probably ones in which the fuel-air ratio is nearly 
stoichiometric. The cooler region near the center (as low 
as ~ 1500 K) is probably fuel rich. Injected fuel has 
probably not penetrated to the sidewalls (~ 1 200 K to 
1300 K). Data are not acquired at planes 6 and 7 in the 
piloted case due to lack of time. However, similarity in the 
pressure distributions between piloted and unpiloled cases 
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Un piloted 


Piloted 


Unpiloted 


Piloted 



unpiloted cases and piloted cases 


downstream of about a = 0.75 m suggest that the 
temperature maps would be at least qualitatively similar, 
except perhaps close to the top wall. 

Consider now the rms temperature surfaces. All the 
surfaces have irregular oscillations of amplitude ~ 100 K 
that are probably simply due to random error. However, 
thete aie signiticant features present. At plane 1 . unpiloted 
case the rms surface varies between ~ 100 K and - 230 K. 
The standard deviations obtained from measurements at 
fixed points varied from 140 K to 220 K. significantly 
greater than that inherent in the instrument 100 K). 
Thus, there are significant real temperature fluctuations in 
the flow exiting the heater, with standard deviation ~ 100 
K or greater. The rms level is comparable in the piloted 
case, except that it rises close to the top surface due to the 
pilot fuel plumes near the top surface. 

Several regions of large temperature rms fluctuation 


levels may be observed. In general, if combustion occurs, 
temperature variations arise due to variation in the local 
equivalence ratio (composition). Temperature variations 
also arise for fixed equivalence ratio due to variation in 
the extent of combustion, from no combustion to complete 
consumption of either fuel or oxidant. Thus, observed 
fluctuations may be due to intermittency in gas 
composition, intermittency in extent of combustion, or 
both. One such region, with peak rms level - 440 K, lies 
between the center of the fuel plume and the top wall, 
plane 5, unpiloted case. A 2 nd such region with peak rms 
level ~ 500 K, lies near the fuel plume, plane 3, piloted 
case. In both cases, these regions appear to be in regions 
where combustion is just beginning or is incomplete, 
suggesting that intermittency in degree of combustion is 
the main cause of the fluctuations. The result for the 
unpiloted case is particularly interesting since there is 
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little indication of combustion in the mean temperature 
plots. One further region of large rms fluctuation level, 
with peak rms level - 525 K, lies on the right side near the 
top wall, plane 5. piloted case. The mean temperature is - 
2200 K at this location (much higher than for the other 
regions discussed), so the level of fluctuations normalized 
by the mean is much lower. It is believed that curve fitting 
accuracy is reduced at the edge of the data plane. Thus, 
this region may simply be explained by random error. 
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Figure 12 Probability density distributions of 
temperature near horizonal centerline, unpiloted. 


Estimates of the probability density functions (pdfs ) 
of temperature for several points near the horizontal 
centerline of each plane are presented in Figs. 12 and 13, 
for unpiloted and piloted cases respectively. These 


AIAA-2002-0743 

estimates are obtained by forming a histogram, with 100 
K wide bins, of each fixed-point data run. The histograms 
are normalized by the number of points in the run and bin 
width. In the figures the vertical axis is temperature and 
each PDF is represented by a vertical stripe centered on 
the v location (horizontal axis) of the point. The color 
indicates probability denoted by the color map. The 
previously discussed regions of high rms levels are visible 
in these figures as regions of broad pdfs. However, it may 
be seen that these pdfs are not always symmetrical. 
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Figure 13 Probability density distributions of 
temperature near horizontal centerline, piloted. 

The fitted mean and rms temperature surfaces, the 
temperature histograms, and the surface-pressure and 
heat-flux data provide a consistent description of the flow r 
field. In the unpiloted case, no significant or only small 
and intermittent combustion of the injected fuel is 
observed ahead of x = 0.5 m. Downstream ot x = 0.75 m 
combustion appears nearly complete. In the piloted case, 
combustion of the pilot fuel appears to take place between 
.v = 0. 1 22 m (at the step) and 0.274 m. There also appears 
to be significant combustion of the main injected fuel by 
this location. These results are not consistent with the 
CFD calculations performed prior to the commencement 
of the experimental work 2 , which predicted combustion in 
the vicinity of injection in both unpiloted and piloted 
cases. Consequently, this experiment provides a test case 
for CFD that is a more challenging than anticipated: 
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accurate calculation will require accurate modeling of the 
chemical kinetics and turbulence -chemistry interactions as 
well as accurate modeling of the turbulent mixing. 

Conclusions 

CARS thermometry and modern design of 
experiments techniques have been successfully applied to 
a simple supersonic combustor. Functions that provide an 
estimate of the mean temperature and mean square 
deviation surfaces at various planes in the flow have been 
obtained. Histograms of temperature are presented. 
Additional data includes surface pressure and 
temperature. Boundary conditions and uncertainties are 
well characterized, and the data are useful for the 
validation of computational fluid dynamics (CFD) codes 
used in the design of supersonic combustors. Accurate 
calculation of the aggregate data set will require that the 
CFD code model both the turbulent mixing, and the 
turbulent kinetics and turbulence-chemistry interactions. 
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